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■ ABSTRACT 

Rapid and strong changes in the Blazhko modulation of RR Lyrae stars, as they 
have recently been detected in high precision satellite data, have become a cru- 
cial topic in finding an explanation of the long-standing mystery of the Blazhko 
, effect. We present here an analysis of the most extreme case detected so far, 

■ the RRab star V445 Lyr (KIC 6186029) which was observed with the Kepler 
' space mission. V445 Lyr shows very strong cycle-to-cycle changes in its Blazhko 

modulation, which are caused both by a secondary long-term modulation period 
as well as irregular variations. In addition to the complex Blazhko modulation, 
\ V445 Lyr also shows a rich spectrum of additional peaks in the frequency range 

. between the fundamental pulsation and the first harmonic. Among those peaks, 

the second radial overtone could be identified, which, combined with a metal- 
licity estimate of [Fe/H]=-2.0 dex from spectroscopy, allowed to constrain the 
mass (0.55-0.65 Mq) and luminosity (40-50 Lq) of V445 Lyr through theoretical 
' Petersen diagrams. A non-radial mode as well as possibly the first overtone are 

\^ , also excited. Furthermore, V445 Lyr shows signs of the period doubling phe- 

nomenon and a long term period change. A detailed Fourier analysis along with 
a study of the 0-C variation of V445 Lyr is presented, and the origin of the 
additional peaks and possible causes of the changes in the Blazhko modulation 
are discussed. The results are then put into context with those of the only other 
star with a variable Blazhko effect for which a long enough set of high precision 
continuous satellite data has been published so far, the CoRoT star 105288363. 
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RR Lyrae stars, which are low-mass helium burning stars 
on the horizontal branch, were long thought to be rather 
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simple radial pulsators. They follow a period-luminosity- 
color relation which makes them valuable distance indi- 
cators, and because of their age and evolutionary status, 
they are also u sed to study t he formation and evolution 
of the Galaxy (|Catelanll2009h . They can oscillate in the 
fundamental radial mode (type RRab), the first overtone 
(type RRc) or both modes simultaneously (type RRd), 
and their high amplitudes of up to 1.5 mag in V for RRab 
type stars made their variability easy to discover, so that 
they have been known since the end of the XIX century. 

Already more than a hundred years ago, however, 
it turned out that there is an a spect to RR Ly rae stars 
which is not understood at all: iBlazhkcJ (|l907l ) found a 
"periodic change in the period" of RW Dra, which he 
could not explain, an d which still remains unexplained 
today. IShaplevI l|l916l ) later found in his observations of 
RR Lyrae that also the brightness of the maxima and 
the light curve shape show periodic changes. With in- 
creasing data quality in the recent past, the unsolved 
problem got even more severe, as it turned out that not 
just a rather small fraction of exceptional RR Lyrae stars 
were affected, but probably around 40-50 per cent of 
all RRab stars (iKolenberg et al.ll2010l : iBenko et al]|2010l : 
ljurcsik et al.l [20091 ). Also among RRc type stars, ampli- 
tude and ph ase modulation w as found to be surprisingly 
widespread (|Arellano Ferro et al.ll201^ V 

This so-called Blazhko effect was long thought to be 
a periodic/regular phenomenon. Traditionally, only one 
Blazhko period was assigned to each modulated star, 
and the phenomenon was expected to repeat in every 
Blazhko cycle, agreeing with the widely used definition 
that "the Blazhko effect is a periodic amplitude and/or 
phase modulation with a period of several 10 to 100 pul- 
sation periods" . There were several reports about changes 
in t he Blazhko modulati o n of v arious stars (see section 
5 of iGuggenberger et al.l l|201ll ) for a recent summary), 
but those reports usually had to rely on sparse data with 
large gaps, so that it was impossible to say when exactly a 
change took place and whether it happened continuously 
or abruptly. The Blazhko effect was therefore still con- 
sidered to be a strictly repetitive phenomenon with only 
some rare exceptions showing secondary modulation pe- 
riods (for example CZ Lac, Sodor ct al. 2011) or changes 
on very long time scales. It was not until the availability 
of ultra-precise data from space missions like CoRoT that 
strong and irregular cycle-to-cycle changes of a Blazhko 
star were documented and that it became obvious that 
seemingly chaotic phenomena need to be accounted for 
when modelling the Blazhko effect. 

While the detection of cycle-to-cycle changes in the 
Blazhko modulation posed a significant challenge for 
all classical models that required a clo ck-work-l i ke be - 
haviour, some new ideas were published. IStotherd l|200d ) 
suggested that transient small-scale magnetic fields mod- 
ulate the turbulent convection inside the helium and hy- 
drogen ionisation zones, a mechanism which certainly 
could explain subsequent Blazhko cycles of differing 
strength. This scenario, however, was re cently tested 
on the basis of hydrodynamical models by ISmolec et al.l 
(|201ll) who found that it was not possible to reproduce 
all observed properties of the light curve, even when al- 
lowing a huge modulation of the mixing length. On the 



other hand, iBuchler fc KolldthI (|201ll ) successfully mod- 
eled both regular and irregular modulations by using the 
amplitude equation formalism. In their models, a strange 
attractor in the dynamics causes chaotic behaviour. The 
9:2 resonance between the fun damental mode and the 
9th overtones that was found bv lKollath et all l|201ll ') to 
be the reason for the recently disco vered phenomenon of 
period doubling l|Szab6 et al.ll2O10l ') , seems to play an im- 
portant role in causing a Blazhko modulation. 

We present here a study of the most extreme case de- 
tected so far: V445 Lyr (KIC 6186029). This RRab star 
not only shows the strongest cycle-to-cycle changes found 
up to now, but also a rich spectrum of additional modes 
in the region between 2 and 4 d"'^. Both Fourier and O-C 
analyses are used to investigate the variability of the pul- 
sation and of the modulation (Section[3|, and the results 
are compared to the case of CoRoT 105288363 in Sec- 
tion [6l Spectroscopy as well as theoretical Petersen dia- 
grams are used to determine the fundamental parameters 
like metallicity, luminosity and mass (Section [S]). Addi- 
tionally, the new analytic modu lation approach for data 
analys is recently proposed by iBenko. Szabo fc Paparol 
l|201ll ) is applied to the data in Section [4.41 



2 BACKGROUND INFORMATION ON 
V445 LYR 

V445 Lyr, with the coordinates RA 18h 58m 26s and Dec 
41° 35' 49" (J2000) is also known as KIC 6186029, or 
GR244, and has a Kepler magnitude of Kp=\7A. Two 
pub lications from t he pre- Kepler era exist for this tar- 
get: [Roman^ (| 19721 ') fou nd it to be variable and classified 
it as an RR Lyrae, and lKukarkin et al.l lll973l) in cluded 
it into his name list of variable stars. iRomand (1972) 
also lists the photographic brightness of maximum and 
minimum to be 15.3 and 17.3 mag, respectively, indicat- 
ing a surprisingly large amplitude of 2 mag, much higher 
than even during extreme Blazhko maxima in the mod- 
ern data. This might at least partly be explained, how- 
ever, by the difference between the observed bandpasses. 
Unfortunately, no details of the observations and no light 
curves are given, and the forthcoming paper that was an- 
nounced by the author could not be found. We therefore 
cannot know if the observed amplitude of 2 mag is real or 
due to possible observational erro rs. No err o r estim ations 
of the observations were given by | Roma no| (|l972l ). 

Since Kepler data have become available, two more 
publications have dealt with V445 Lyr, both p resent- 
ing the Kepler data up to Q2: ISzabo et al.l (|2010l ) listed 
V445 Lyr as a possib le candidate for th e period dou- 
bling phenomenon, and lBenko et al.l (|2010l l already noted 
changes in the Blazhko modulation of V445 Lyr, and re- 
ported the presence of radial overtones. 



3 KEPLER PHOTOMETRY 

The Kepler space mission was launched o n March 6, 
2009 into an Earth-trailing heliocentric orbit l|Koch et al.l 
I2OIOI ). Its primary purpose is the detection of Earth-sized 
planets in the habitable zone of solar-like stars through 
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Figure 1. Panel a shows the complete light curve of V445 Lyr presented in this paper, including quarters Q1-Q7. Due to failure 
of a module, there are no data from Q5, resulting in the gap from MJD 55273 to 55372. Panel b illustrates the extreme Blazhko 
modulation by comparing a pulsation cycle at Blazhko maximum to one at Blazhko minimum. The light variation has an extremely 
low amplitude during this Blazhko minimum and shows a double maximum. In panel c, a zoom into a region with period doubling 
is shown as an example, and panel d emphazises the distorted light curve shape that occurs during some Blazhko minima. 



the transit method, which requires continuous and ultra- 
precise photometry of over 150000 stars for at least 3.5 
years. This is also the duration of the primary mission. 
Kepler therefore provides not only the longest continuous 
data sets ever observed for RR Lyrae stars, but also does 
so with the highest photometric precision ever obtained, 
as a consequence greatly improving our knowledge about 
stellar pulsations. 

The Kepler spacecraft carries a Schmidt telescope 
with an aperture of 0.95 m and 42 science CCDs 
which cover a field of view of about 115 square degrees 
(| Jenkins et al.ll2010l '). The photometric bandpass ranges 
from 423 to 897 nm, thus avoiding the Call H and K lines 
in the blue, and fringing due to internal reflection in the 
red. The Kepler band is therefore slightly broader than 
a combination of Johnson V and R, and Kepler magni- 
tudes are usually about 0.1 mag from R for most stars 
(|Koch et al.ll201Cl ). 

Every quarter orbit, the spacecraft is rotated in order 
to keep the solar panels oriented towards the Sun, and 
the ra diator wh ich cools the focal plane towards deep 
space (H aas et al 2 01Q ). Data from different quarters are 
denominated Ql, Q2, etc. Each of the quarters which are 
used here has a time base of about 90 d, except Ql which 
covers about 33 d. 

Each measurement is based on a 6.02 s exposure plus 



a 0.52 s readout time. To obtain Long Cadence data (LC, 
29.4 min), 270 measurements are coadded; for Short Ca- 
dence data (SC, about 1 min), 9 exposures are coadded 
^Christiansen et al.llioTll '). 

The time of each measurement is given in truncated 
barycentric Julian Date (HJD-2400000), and refers to 
the midpoint of the measurement. 



3.1 The V445 Lyr data set 

Figure [T] illustrates the data obtained for V445 Lyr. In 
this paper, we present LC data obtained in the quarters 
1 to 7, with a total time base of 588 d. Because of the loss 
of module 3, which happenend in January 2010, there are 
no data available for this target in Q5. Some smaller gaps 
are also present in the data as can be seen in Fig. [T] Some 
of them are due to unplanned safe mode events or loss of 
fine point control, others are caused by the regular down- 
links where the spacecraft's antenna is pointed towards 
Earth for data transmission, and science data collection 
is interrupted. 

Kepler provides light curves in two different formats: 
raw fiux or corrected fiux. The latter has been processed 
for planet transit search by the PDC (Pre-search Data 
Conditioning) pipeline which is known to sometimes re- 



© 2002 RAS, MNRAS 000.[TlfT9l 



4 E. Guggenberger, K. Kolenberg, J. M. Nemec, R. Smolec et al. 



Table 1 . The scaled and detrended data set of V445 Lyr that 
was used for the analysis in this paper. Column 1 gives the 
truncated barycentric Julian date, column 2 the magnitude 
with the average shifted to zero, and column 3 gives the quar- 
ter in which data were obtained. The full table is available 
online only. 



HJD-2400000 


mag 


quarter 


54964.51211 


0.044188798 


Ql 


54964.53254 


0.038912686 


Ql 


54964.55298 


0.017979982 


Ql 


54964.57341 


-0.023960246 


Ql 


54964.59385 


-0.099222117 


Ql 



move astrophysical features and which does not preserve 
all stellar variability. Hence here, and for the study of 
variable stars in general, we make use only of the raw 
time series. 

As the spacecraft is rotated every quarter orbit, the 
target falls on different CCDs after each "roll", resulting 
in differences in average flux due to different sensitiv- 
ity levels. Additionally, trends might occur due to image 
motion on the CCD or sensitivity changes. The scaling 
and detrending which has to be performed before start- 
ing the analysis is a delicate task especially for targets 
like RR Lyrae stars which have high ampli tudes, long 
periods and changes in the amplitudes (see ICelik et al.l 
for a detailed discussion). Here, we removed lin- 
ear trends which were determined from a running average 
separately for every quarter, and scaled every quarter to 
the same mean brightness. The continuity of the upper 
and lower envelope of the lightcurve was a good indicator 
of correct scaling. Note that due to the gap during Q5, 
the continuity could not be checked there. For this rea- 
son, the scaling of the data obtained after the gap might 
not exactly be consistent with the scaling applied before 
Q5. 

The scaled and detrended data that were used in this 
analysis are available as online material in the format 
displayed in Table [T] 

Panel a of Fig. [T] shows the complete combined data 
set of V445 Lyr. A total of 5.5 and 3.5 Blazhko cycles 
were observed before and after the large gap, respec- 
tively, and 900 pulsation cycles were observed. Already 
from looking at panel a it is obvious that the Blazhko 
effect shows strong and fast changes, with almost all ob- 
served cycles having a different appearance. Panel b com- 
pares in a phase diagram directly a light curve from a 
Blazhko ma^dmum with one obtained during the subse- 
quent Blazhko minimum, illustrating the extremely low 
pulsation amplitude that occurs during some, but not all, 
modulation minima. Note that due to the very small am- 
plitude and the distorted light curve showing a double 
maximum (see also panel d of Fig. [1]) it would not be 
possible to recognize the RR Lyrae nature of this star 
when observing it only during Blazhko minimum where 
the peak-to-peak amplitude can be as low as 70 mmag. 
Extremely low pulsation amplitud es during Bla zhko min- 
ima were recently also reported bv lSodor et al.l (2012) for 
two of RRab type stars, and the occurrence of a strong 
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Figure 2. Spectral window function of the complete V445 Lyr 
data set. Lower panel: magnification of the region between 
5 and 7 d~^, showing some examples of the low-amplitude 
features caused by the reaction wheel. 



bump during Blazh ko minimum was noted in RZ Lyr 
(|jurcsik et al.ll2012l ). In panel c of FigH] a small part of 
the light curve around JD 2455386 is shown, illustrating 
the phenomenon of period doubling (alternating higher 
and lower light maxima) which was recen tly discovered 
in Kepler Blazhko stars (jSzabo et al.ll2O10h and which is 
also present in V445 Lyr. 



4 ANALYSIS AND RESULTS 

Due to the quasi continuous coverage and the high photo- 
metric precision, the conditions for a Fourier analysis are 
more than favorable. The spectral window function (see 
Fig. ^ is almost perfect, without any alias peaks visible 
at first glance. Because of the Earth-trailing orbit, no or- 
bital frequencies like the ones that cause trouble in the 
data from many other space missions are present here. 

A well-known feature in the Kepler data are the 
momentum desaturations of the reaction wheel, which 
happen every 2.98 d. During the thruster firings which 
are necessary to release the angular momentum that has 
built up in the reaction whe els, the space craft momentar - 
ily loses fine point control (|Van Cleve fc Caldwelll [20091 ). 
This results in a missing data point during every desatu- 
ration in the otherwise regularly spaced data, leading to a 
comb-like structure in the Fourier transform with a spac- 
ing of 0.335 d~^ . In the spectral window, those peaks have 
a very low normalized amplitude of only 0.004 mag. A 
zoom into the spectral window (see lower panel of Fig. [S]) 
reveals this comb of tiny peaks. When carefully remov- 
ing all high-amplitude peaks from the Fourier spectrum 
before turning to interpreting features with amplitudes 
which are orders of magnitude smaller as are discussed 
in section 14.1.51 those alias peaks are not expected to 
cause any trouble. 

Other features in the spectral window function are 
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two harmonics of the monthly data downlink frequency, 
at 0.065 and 0.13 d~^, with normalized amplitudes of 
0.05 and 0.03 mag, respectively. 

Due to the time base of 588 d, the Rayleigh frequency 
resolution is 0.0017 d~^. 



4.1 Fourier analysis of the full data set 

As a first step, a Fourier analysis of the complete data 
set was performed, keeping in mind that changes of the 
fundamental period might occur during the time of the 
observations, and that the changes of the Blazhko modu- 
lation cause a large number of peaks close to the classical 
pulsation and modulation components. The Fourier anal- 
ysis of the full data set, however, is necessary to obtain an 
overall picture of the pulsation and modulation proper- 
ties of V445 Lyr, to find the mean values of the pulsation 
and modulation periods and to detect possible long-term 
periodicities which can be resolved with a long time base 
only. 

The Fourier analysis of the complete dat a set was 
performed with Period 04 (iLenz fc Bregerll2005l ) and then 
checked with SigSpec (|Reegen 2007 ). The results agreed 
within the errors. The Fourier analysis revealed a mean 
pulsation period of 0.513075 ± 0.000005 d and a mean 
Blazhko period of 53.1 ± 1 d with the ephemeris 



HJD (Tmax, pulsation 

) = 2455550.514 + 0.513075 x E 

pulsation, 

HJD (Tn,ax,Blazhko) = 2455534.2 + 53.1 X Esiazhko. 

Note that the presence of close peaks has a great influence 
on the fitting and frequency optimization procedure, and 
the errors are therefore larger than would normally be 
expected for a Fourier fit to data of the given quality. 

4.1.1 Multiplet components 

The Blazhko multiplets (i.e., the pattern of peaks which 
is typical for Blazhko stars with peaks at the positions 
kfo ± rifs, with k and n being integers denoting the 
harmonic order and the multiplet order, respectively, 
and with /o and /s denoting the fundamental and the 
Blazhko frequency) were found to be very asymmetric 
in amplitude. Much higher amplitudes appeared on the 
right (higher frequency) side than on the left. This is the 
more common case, which is ob served in thr ee fourths 
of all Blazhko RRab stars (Alco ck et al.|[200^ '). Compo- 
nents were detected up to quintuplet order on the right 
side, while on the left side of the main pulsation compo- 
nent, only one side peak (i.e., a triplet component, n = 1) 
could be found. In some orders, peaks were found near 
the positions that would be expected for septuplets, but 
the deviations from the exact frequency values were con- 
sidered too big to identify the peaks safely as septuplet 
components. 

In addition to the classical Blazhko multiplets 
around the fundamental mode and its harmonics, also 
a secondary modulation component, hereafter fs, could 
be identified. It manifests itself as a series of additional 
peaks next to the classical modulation sidepeaks of the 
Fourier spectrum, appearing in every order with a spacing 



of 0.00698±0. 00027 d~^, indicating a secondary modula- 
tion period of 143.3 ±5.8 d. The secondary peaks appear- 
ing on the right side of the classical peaks show surpris- 
ingly high amplitudes, while the secondary peaks to the 
left of the classical multiplet have small amplitudes and 
could only be detected after prewhitening a large num- 
ber of higher amplitude peaks. Interestingly, the peaks 
belonging to this additional series could also be detected 
in higher multiplet orders (n = 3) than the classical mul- 
tiplet, making it difficult to explain the peaks as com- 
bination frequencies. Their amplitudes seem to decrease 
less quickly with increasing harmonic order k (see also 
section |4.1. 3 |l . making them easier to spot at high orders. 

Figure [3] shows the Fourier transform of the data, 
providing also zooms into the regions around the 3'^'*, 4*'' 
and 5**^ harmonic order where both the classical mul- 
tiplet and the additional components can be seen. The 
highest peaks among the multiplet components, which 
can be seen even before prewhitening the original data, 
are marked with arrows. 

Figure 13] illustrates the pattern of detected frequen- 
cies in the vicinity of the fundamental pulsation and 
its harmonics in the style of an "echelle" diagram. In 
this diagram, which is similar to the diagrams used to 
unveil equally spaced peaks in helioseismology, the fre- 
quency of each peak is plotted against (/ modulo /o), 
i.e., f/fo - INT{f/fo), or f/fo - INT{f/fo) - 1 for 
peaks to the left of the harmonic, therefore clearly re- 
vealing patterns which repeat in every harmonic order. 
Peaks belonging to the same group of combinations align 
in vertical ridges. We stress that unlike in the helioseis- 
mic application, where the ridges denote different radial 
orders of same degree, in this case the echelle diagram 
only serves the purpose of displaying in a very practical 
and easy way the repeating patterns in different harmonic 
orders of non-sinusoidal fundamental radial pulsation. 



4.1.2 Devtatton of the harmonics 

Due to the non-sinusoidal light curve shape typical for 
RR Lyrae stars, harmonics of the fundamental mode 
are expected to appear at the frequencies kfo, where k 
is an integer denoting the harmonic order. The classi- 
cal Blazhko multiplets in modulated stars are spaced 
equidistantly, implying frequency values of kfo ± nfs 
with n denoting the multiplet order. Long-term period 
changes and close peaks caused by irregular phenomena, 
however, can distort this frequency pattern. When 
analyzing time series of Blazhko RR Lyrae stars, there 
are two options for fitting the data: one is to fix the 
frequencies of the harmonics and Blazhko multiplets to 
their expected values of kfo + nfB, therefore reducing the 
number of free parameters in the fit. The other option 
is to let all parameters, including the frequency values, 
free. When the latter option was applied to this data set, 
the harmonics were observed to deviate systematically 
and significantly from their expected values, which can 
also be noticed as a slight rightwards tilt of the ridges 
in the echelle diagrams (Fig. [4]). Normally, one would 
expect that this is simply caused by a wrong value for 
fo, but in this case, no value of fo could be found which 
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Figure 3. Fourier transform of the complete V445 Lyr data 
set (upper panel). In the lower panels (b, c and d), the regions 
around the harmonics fc/o with k=3, 4 and 5 are shown in de- 
tail. Black arrows indicate the peaks belonging to the classical 
Blazhko multiplet with fc/o — /s, fc/o, kfo + fs and fc/o + 2/B, 
while grey arrows indicate some of the peaks belonging to the 
secondary multiplet which within itself has a spacing equal 
to the Blazhko frequency, but is shifted with respect to the 
classical multiplet by 0.00698 d~^. Note that from order k=4 
onwards, the highest amplitude peak in this order is a peak 
belonging to the secondary multiplet 
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Figure 4. "Echelle" diagram of the peaks detected in the 
vicinity of the fundamental pulsation mode and its harmon- 
ics. Frequency is plotted versus / modulo /o, so that regu- 
larly spaced peaks which appear in every order align in verti- 
cal patterns, making it easy to identify combinations with /o- 
Peaks belonging to the category fc/o it nfg are shown as black 
squares, peaks of the group fc/o it nfg + fg are shown as grey 
circles, components of fc/o -t nfB + 2fs are plotted with light 
grey triangles. Open circles denote unidentified and/or unre- 
solved peaks, some of which originate from the non-repetitive 
nature of the modulation and long-period changes. 



could solve the issue, i.e., every detected harmonic, when 
divided by its order, required a difTerent fo- We therefore 
decided in favour of the more pure approach and did not 
fix the frequency values to the expected positions, but 
left all parameters free in the fit of the complete data 
set. The problem disappeared, however, when analyzing 
different subsets of data separately (see section [4^ . and 
we therefore suspect it to be either the result of period 
changes which take place during certain seasons (see 
also section 14. 3p and/or close unresolved peaks which 
are known to strongly infiuence the results of both the 
Fourier analysis and the multisine fitting procedure. 



4-1.3 Amplitudes versus harmonic order 

It is a well-known fact that in Blazhko RRab stars the 
amplitudes of the multiplet side peaks decrease less 
rapidly with harmonic order than t hat of the ma i n com - 
ponent. This was first described bv lJurcsik eFal] (|2005l ). 
and th en confirmed for other well-studied stars like 
SS Fo r (|Kolenberg et al. 2009 ). RR Lyr (.Kolenberg et all 



„ _ _ ^ , _ rg 

I2OIII ) and CoRoT 105288363 l|Guggenberger et al.ll20lil ) 
As expected, the result is the same for V445 Lyr (see 
Fig. [S]), but in an extreme way with the ampitude 
of the right side peak exceeding that of the main 
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Figure 5. Amplitudes of all components versus harmonic or- 
der. Components belonging to the classical Blazhko multiplet 
are shown with green symbols, peaks which are combinations 
with fs are shown in orange, and components with a 2/s term 
are shown in magenta. It can be seen that the amplitude of 
the fs multiplets decrease less steeply than those of the clas- 
sical Blazhko multiplet, reaching the same amplitude in the 
4th order and dominating the frequency spectrum for higher 
orders. 



component as early as in the 3rd order. Additionally, 
the amplitudes of the secondary modulation multiplet, 
i.e. combinations with the secondary modulation fs, 
could be studied in this case. It turned out that the 
amplitudes of the secondary multiplet components 
decrease even less steeply, therefore dominating the 
Fourier spectrum from the 4th order onwards. The 
strongest signal then comes from the combination 
fo + /b + fs, i-e. the peak on the right side of the right 
triplet component. Also, combinations with the term 2/s 
could be detected, and their amplitudes are also plotted 
in Fig.O They also show a very slow amplitude decrease. 



4-1-4 Number of relevant frequencies 

Due to the dense spectrum of peaks which is caused by 
the cycle-to-cycle changes of the Blazhko effect, an anal- 
ysis in the classical sense, i.e. taking into account all fre- 
quencies down to a certain S/N level or a certain sig- 
nificance criterion, might not be optimal in a case like 
V445 Lyrae, as it does not yield meaningful results. A 
large number of the detected peaks is likely to be the 
result of "stellar noise", caused by irregular and/or long- 
periodic phenomena, and many of them are not resolved 
with the available time span. Tests revealed that as many 
as 771 frequencies can be found when performing an anal- 
ysis until the generally adopted criteria of significance are 



reached. Many of them were not resolved, and many could 
not be attributed to any combination of other modes, 
and did not show repeating patterns in the echelle di- 
agrams. Therefore, instead of choosing the classical ap- 
proach, the analysis was stopped after a certain number 
of the highest peaks in every harmonic order was found 
and prewhitened, usually around 20 peaks per order. It 
turned out that after subtracting approximately 20 peaks 
in a given harmonic order, no meaningful combinations 
could be identified among the following peaks, and many 
unresolved peaks appeared. In Fig. [3] only the highest 
peaks of every order are shown, already including some 
unresolved peaks which could not be avoided due to their 
high amplitude. As the Nyquist frequency of LC data is 
24.4 d~^, 12 harmonic orders could be observed, and 239 
frequencies were subtracted around the main pulsation 
components until the attention was turned towards the 
additional peaks which are present in the region between 
the harmonics of the fundamental mode (see next sec- 
tion). 

Altogether, 239 frequencies were included, of which 
one is the fundamental mode, 9 are harmonics of the 
fundamental mode, 104 are combinations of fo or its 
harmonics with the Blazhko frequency fs and/or the 
secondary modulation fs, and 125 are unidentified. 



4-1-5 Additional frequencies: radial overtone and 
non-radial pulsation 

It was already noted by iBenko et al.l (|2010l ) that 
V445 Lyr shows a rich spectrum of frequencies in the re- 
gion between the harmonics of the fundamental mode. 
Those frequencies are not at all typical for ab-type 
RR Lyrae stars and have never been detected in such a 
large number in an RR Lyrae st ar. Some peaks wer e sus- 
pected to be radial overtones bv lBenko et al.l l|2010l ), and 
also frequencies which are half-integer multiples of the 
fundamental mode are expected to appear in this region 



as a consequ ence of the period- doubling phenomenon as 
described bv lSzabo et al.l l|2010l ). But there is more than 
this to V445 Lyr. 

Fig. [6] shows the frequency spectrum after subtrac- 
tion of the relevant peaks around the multiples of the 
fundamental pulsation as discussed in the previous sec- 
tion (their former positions are marked with arrows) . The 
additional frequencies can clearly be seen to be the domi- 
nant signal with an amplitude of 3.7 mmag for the highest 
peak. The region between the fundamental mode and its 
first harmonic is indicated with a grey box and enlarged 
in the insert (panel b). Four frequency regions with en- 
hanced signal can be noted in the enlargement: around 
2.65, 2.8, 2.9 and 3.33 d~^. This pattern repeats in ev- 
ery harmonic order, indicating combinations of the fre- 
quencies with the fundamental mode and its surrounding 
peaks. The presence of combinations is a strong evidence 
that the additional signal is not introduced by a possible 
background star. 

In a Fourier analysis of the relevant frequency re- 
gions, 80 peaks (including combinations) were consid- 
ered siginficant and subtracted. A closer inspection of 
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Figure 6. Fourier spectrum after prewhitening 239 frequen- 
cies in the vicinity of the fundamental mode and its harmon- 
ics. Very clearly, additional frequencies can be seen between 
the harmonic orders, as well as their combinations with /o 
which are present up to high harmonic orders. Arrows mark 
the places where the harmonics and the Blazhko as well as 
the secondary multiplets were located before prewhitening. 
Some signal remains around their positions, as discussed in 
section [4. 1.41 The insert (panel b) shows a zoom into the re- 
gion between the fundamental mode and the first harmonic 
(2-4 d-i), indicated as a grey square. 



the result revealed that the dominant peaks formed com- 
binations not only with the fundamental mode but also 
with the Blazhko multiplet peaks (including quintuplets!) 
and in some cases also with the peaks belonging to the 
secondary multiplet. Also negative combinations such as 
In — fo — fs occur. Significant combination peaks can 
be traced up to the fifth harmonic order, but an excess 
in signal is visible in the Fourier spectrum even at much 
higher orders (see Fig. [6]). In Fig.[7l an echelle diagram is 
plotted for the additional peaks, clearly showing the com- 
binations with the fundamental mode aligned in vertical 
patterns. Shaded boxes indicate the typical regions in 
which overtone modes and half-integer combination fre- 
quencies (HIFs) would be e xpected to be situat ed. Please 
note that it was shown by ISzabo et al.l (|2010l ) that due 
to the on- and offset of the period doubling phenomenon, 
the HIF peaks are not necessarily located at the exact 
positions of the half-integer multiples, but might devi- 
ate by several per cent. Therefore, the shaded box at (/ 
modulo /o)=0.5 in Figure[7]has a distinct width. 

In V445 Lyr, the frequency at 2.9256 d~^ deviates 
by only 0.0021 d~^ (i-e., 0.07 per cent) from the exact 
value of 3/o/2. Given the fact that clear signs of period 
doubling are indeed visible in the light curve (see panel c 
of Fig. [T]), and considering the above-mentioned findings 
of ISzab6 et al] l|2010t ). it is quite safe to interpret this 
frequency as an HIF caused by period doubling. We 
hereafter refer to it as fn- Four significant HIFs (which 
can also be interpreted as combinations of fn with /o) 
were found in the Fourier spectrum: 3/o/2, 5/o/2, 9/o/2, 
and 11/0/2. Combinations with the Blazhko frequency, 
both positive and negative, could also be identified (see 
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Figure 7. Echelle diagram of the additional frequencies, show- 
ing how the patterns repeat in every harmonic order. Shaded 
boxes indicate the ranges which would be typical for the first 
and second overtone (Ol and 02, respectively), as well as for 
the half-integer frequencies (HIFs). 



Table [2] for a complete list). 

Another interesting feature is the peak at 3.3307 d~^ 
which shows a frequency ratio of /0//2 =0.585 with the 
fundamental mode, and which we hereafter refer to as 
/2. Its period ratio is typical for the second overtone. 
Peaks with similar period ratios have already been re- 
port ed for several RR Ly rae stars. The first to find them 
were IPoretti et al.l (HoTo) in the CoRoT star 101128793. 
The frequency ratio in their study was 0.584 with the 
fundamental mode at /o = 2.11895 d~^. They interpreted 
the peak as the second radial overtone. The same authors 
also reanalysed the data of V1127 Aq l ijChadid et al.l 
I2OIOI ) and MW Lyr (|jurcsik et al.ll2008[ ) and found fre- 
quency ratios of 0.582 (/o = 2.8090 d"^) and 0.588 
(/o = 2.5146 d~^), respec t ively. In the sample of Kepler 
RRab stars, iBenko et al.l l|2010l ) reported the presence 
of the second overtone in four different stars: V354 Lyr 
(a Blazhko star with /o = 1.78037 d"^), V2178 Cyg (a 
Blazhko star with /o = 2.05423), V445 Lyr (the sub- 
ject of this paper) and the non-modulated RRab star 
V350 Lyr (/o = 1.68282) which was the first example 
of a non-Blazhko double mode RR Lyrae star with the 
fundamental (F) and the s econd overtone (02) excited. 
ICuggenberger et al.l ()201lh found evidence for the sec- 
ond overtone in CoRoT 105288363, a Blaz hko star with 
rapid changes in the Blazhko effect, and iNemec et al.l 
(I2OII') found the second overtone in KIC 7021124, there- 
fore providing another example of a non-modulated RRab 
star pulsating in both F and 02 (/o = 1.606445 d"-^, 
/0//2 = 0.593). 

It is interesting to note how different the stars are 
for which the second overtone has been documented so 
far: Their fundamental frequencies range from about 1.6 
to 2.8 d~^, covering almost the full bandwidth of RRab 
pulsation, and with respect to stability they range from 
non-modulated stars with almost perfectly regular RRab 
pulsation (V350 Lyr, KIC 7021124) to Blazhko stars 
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Lyrae stars in the Galactic bulge ha ve /p //i values 
going down to about 0.726 only (see ISoszvnski et alj 
|201l|). Extremely high metallicity values would be 
necessary according to models ( Smolec fc MoskalikI 
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Figure 8. Amplitudes of the additional frequencies and their 
combinations with /o versus frequency. While the amplitudes 
of and /i decrease exponentially with harmonic order, the 
amplitude of /2 remains almost stable. Amplitude errors are 
smaller than the symbols. 



with a rather regular Blazhko effect (CoRoT 101128793 
and V1127 Aql) and finally to modulated stars that 
show dramatic changes of their Blazhko modulation 
(CoRoT 105288363 and V445 Lyr). Also, they cover a 
significant range of Blazhko periods, from 16.6 d to more 
than 200 d, as estimated for V2178 Cyg. 

The combinations of /2 in V445 Lyr deserve some 
special attention. While as many as 32 peaks are 
detected near the positions of (2 + kfo, and clearly an 
excess of signal is visible in every harmonic order (see 
Figs. [S] and [T]), it was not possible to identify most of the 
detected peaks as exact combinations with the known 
frequencies. For fn (which was discussed in the previous 
paragraph), 12 out of 16 peaks could be attributed to 
combinations while for /2, only 3 combinations with 
fo were found at their exact positions. All the other 
peaks in the vicinity of the combinations deviated too 
much from the calculated values to be safely matched 
with combinations. This is especially remarkable as 
the amplitudes of those peaks are surprisingly large in 
higher harmonic orders compared to the combinations 
of the other additional frequencies. From Fig. [6] it is 
obvious that in the first harmonic order, /2 has a small 
amplitude compared to the other additional frequencies, 
while at the orders 4-5 they become equal, and at higher 
orders, the peaks in the area around /2 + kfo are the 
dominant feature. This is also illustrated in Fig. [S] where 
the amplitudes of the safely identified peaks are shown 
versus frequency. The large number of unidentified peaks 
around /2 might indicate that the amplitude of /2 is 
variable, either irregularly or on a timescale other than 
the Blazhko frequency. This will be discussed in more 
detail in Section [4.2.11 

The highest-amplitude peak among the additional 
frequencies (3.7 mmag) is the one at 2.7719 d^^. 
This peak was in terpreted as the first overtone by 
iBenko et al] (|2010l '). but its ratio with the fundamental 
(/o//jv=0.703) is very low compared to the canoni- 
cal value of 0.74-0.75. Note that the OGLE III RR 



2008 : iPopielski. Dziembowski fc Cassisil I2OO0I : 
iSzabo, Kollath fc BuchleJ |2004| ) to fit this frequency 
ratio with the first overtone. New spectroscopic results 
revealed, however, that the metallicity of V445 Lyr is 
[Fe/H]=-2.0 ± 0.3 (see Section[5]), rendering it impossible 
to explain this frequency with a radial overtone mode. 
/at is there fore most likely to be a non-radial mode. We 
note that IVan Hoolst. Dziembowski fc KawaleJ (| 19981 ) 
found in their non-adiabatic non-radial calculations 
the excitation of non-radial modes in the vicinity of 
the radial mode in RR Lyrae st ars to be very likely, 
and [Dziembowski fc Cassisil (| 19991 ) noted in their model 
survey the presence of strongly trapped non-radial modes 
with very high growth rates near the first overtone. 

Among the additional modes in V445 Lyr, /at is 
the one that shows the clearest pattern of combination 
frequencies: of 20 peaks which were found significant in 
relation to /at, all 20 could be unambiguously identified 
as exact combinations with fo, the Blazhko frequency 
and the secondary modulation frequency (see also 
Table [5] and Fig. [S] which illustrates the regular pattern 
of combination peaks). We note that the possible n on- 
radial mode which was found bv lChadid et al.l (|2010h in 
VI 127 Aql has a very similar frequency ratio (0.696), 
maybe hinting at a possible systematic preference in 
non-radial mode excitation in RR Lyrae stars. 

There remains, however, the fourth region of in- 
creased signal with a main peak at 2.6676 d^^, which, 
with a frequency ratio of /o//i=0.731 is in principle in 
the possible range of the first overtone pulsation. One 
has to note that double mode RR Lyrae stars usually 
follow a well-defined empirical sequence in the Petersen 
diagram, in oth er words, there is a relation between /p 
and /0//1 (see Popielski. Dziembowski fc Cassisil I2OO0I : 



ISoszvnski et aLlbOOgl ). If the peak at 2.6676 d~^ is indeed 
the first overtone, V445 Lyr would be an exception to this 
relation, which is very unlikely. On the other hand, out- 
liers from the se quence similar to V445 Lyr have recently 
been reported bv lSoszvnski et all (|201ll ') in the OGLE III 
survey of the Galactic Bulge. 

The metallicity needed to reproduce a frequency ra- 
tio of /o//i=0.730 with models (Z=0.004, see Figure [TO)) . 
is much larger than the spectroscopic value (Z=0.0002, 
see also Section [SJ. Also, the average frequency values 
are quite far (2/i — fo ~ f2 ~ 0.055) from the resonance 
condition that could explain the presence of the second 
overtone by resonant excitation, and which would also 
have the power to shift the frequency away from the 
expected value in the Petersen diagram. We note, 
however, that the frequency values of /i and /2 are not 
strictly constant during the observed time span, but 
undergo irregular fluctuations. We performed a time- 
dependent Fourier analysis (see also section 14.2. ip and 
found the resonance criterion to be fulfilled occasionally. 
For /i, one combination with fo and 5 combination 
frequencies with fo and fs were found, leaving the other 
6 significant peaks unidentified. One of them, a peak 
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Table 2. List of the highest peaks in every one of the four 
regions of excess signal, and their combinations. 
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Figure 9. Detail of the echelle diagram, showing only the 
vicinity of the mode /n. All detected peaks could clearly 
be identified as combinations of with either the classical 
Blazhko multiplet or the secondary multiplet. The fact that 
no unidentifiable peaks are among the highest ones points to- 
wards a very stable amplitude of this mode, compared to the 
other additional peaks that were found in this star. 
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at 2.639 d^^ (with a frequency ratio of 0.739 with the 
fundamental mode), would fulfill the resonance criterion, 
but its amplitude is only 1.9 mmag, compared to 2.8 
mmag of fi . We therefore conclude that the identities 
of the peaks at 2.6676 d~^ and 2.639 d~^ cannot be 
unambiguously assessed. 

The highest peaks in every region, as well as their 
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Figure 10. Petersen diagram for the peak at 2.6676 il- 
lustrating that a metallicity of Z=0.004 would be necessary to 
identify this peak as the first radial overtone mode. Models 
around the spectroscopically derived value of Z=0.0003 (see 
Section [SJ are also shown for comparison. Different symbols 
indicate different luminosities, and different line styles indicate 
different masses (see also Section O . Models were calculated 
with the Warsaw codes ISmolec fc Moskalik.i2008 ). 



combinations, are listed in Table (2] Altogether, 80 fre- 
quencies were found to be significant in the vicinity of the 
additional peaks: four of these were independent frequen- 
cies, 40 were combinations of these independent terms 
with fo, /fl and/or fs, and 36 were peaks which could 
not be identified as combinations. Together with the 
239 frequencies found in the vicinity of the fundamental 
mode and its harmonics (see Section I4.1.4p . this re- 
sults in a total of 319 frequencies included in the analysis. 
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4.2 Separate analysis of subsets 

To study the time-dependent behaviour of the frequency 
pattern, we used bot h Period04 an d the time-resolved 
feature of Sie Spec (|Reegenl |2007| ). which makes it 
possible to analyse large numbers of data subsets in an 
automated way. The top panel of Fig. [TT] illustrates the 
variation of /o with time as calculated with SigSpec 
for overlapping subsets of 15 d. With Period04, we 
analyzed larger subsets of data, consisting of one or more 
quarters each. During this process we noticed that in the 
subset containing Ql to Q4 (Blazhko period 53.9 d), no 
deviation of the harmonics as described in section 14.1.21 
is observed. This phenomenon seems to occur only in 
quarters Q6 and Q7 were the Blazhko periods found in 
separate analyses of the quarters were 79.8 and 80.4 d, 
respectively. 



4.2.1 Stability of the additional frequencies 

The results discussed in Section [4. 1 . 5 1 hint towards a vari- 
ability of the amplitude of /2, and we therefore studied 
the temporal evolution of this peak in detail. Using the 
time-resolved mode of Sigspec, we performed a Fourier 
analysis of overlapping bins with a duration of 15 d in 
steps of 2 d, limiting the frequency range to the region 
around /2. The resulting amplitudes are plotted in panel 
b of Fig. 111! While there is a clear variation of the ampli- 
tude ranging from 2 mmag to about 7 mmag, no clear pe- 
riodicity is discernible. We also performed a Fourier anal- 
ysis on the resulting amplitude curve and found no signif- 
icant frequency of variability. We note that in Fig. 1111 all 
values of A{f2) were included, regardless of significance 
of the frequency. Only bins with an insufficient number 
of data points and/or bad frequency resolution were dis- 
carded. Therefore we provide as a supporting plot the 
time-dependent significance of the peaks in panel c, where 
the most commonly used significance criterion of sig^5 
is indicated with a dashed line. 

For better orientation, the bottom panel shows the 
full data set (light curve) of V445 Lyr, and shaded boxes 
indicate the regions of enhanced /2 amplitude. There 
seems to be a preference for phases close to Blazhko 
minima, which are also the phases were the minima of 
the funadamental frequency /o occur, but there is no 
strict rule that can predict a high /2 amplitude. This 
irregular amplitude variability can explain the numerous 
peaks around /2 and its combinations which can be seen 
in the echelle diagram (Fig. [7]) and which were discussed 
in section 14.1.51 

We also checked the temporal variations of the 
frequency values of all the additional frequencies and 
found, in addition to the variation of /o which is plotted 
in panel a of Fig. 1111 slight irregular variations of /i 
and /2, which lead to an exact parametric resonance 
during some time intervals in the observed data. This 
might result in the momentary and transient excitation 
of /2 which is seen in Fig. 1111 We note that from the 
theoretical point of view, a resonance is necessary to 
excite the second overtone in this parameter range, be- 
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Figure 11. Variation of /o over the complete data set (panel 

a) . Middle panels show the variation of the amplitude (panel 

b) and the significance (panel c) of the second overtone /2, 
calculated for bins of 15 d duration with a step width of 2 d. 
The bottom panel (d) shows the light curve for orientation. 
The times of significant /2 amplitude are marked with shaded 
boxes. 



cause otherwise this overtone would not become unstable. 



4.2.2 Variation of the Blazhko modulation parameters 

The modulation parameters are normally used to de- 
scribe the properties of the Blazhko modulation of a given 
star. The traditional parameters are Rk = A+/A- where 
j4+ and A- are the amplitudes of the peak on the higher 
frequency and lower frequency side of the triplet, respec- 
tively, and the phase difference /S.ipk = V'+ ~ The 
parameter k denotes the harmonic order. Also the asym- 
metry parame ter Q = (A+ — A- ) / {Aj^ -\- A-), which was 
introduced bv lAlcock et~al] (j2003l ). is widely used, as well 
as the power difference = AX — A?_ . Th i s par am- 

eter was recently shown by ISzeidl &: Jurcsilj (|2009l ) to 
be the physically most meaningful one, as it is directly 
correlated to the phase difference between the amplitude 
and the phase modulation components in their model of 
modulated oscillation. 

As the main aspect of V445 Lyr is the variability of 
the Blazhko cycles, we show here not only the average 
Fourier parameters for every order (which are listed in 
Table [3]) but also their time dependent behaviour in Fig- 
ure [121 To calculate the modulation parameters for every 
cycle, the data set was divided into bins, starting and 
ending around Blazhko minimum. The length of the bins 
was about 60 d, with the exception of the last bin which 
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Figure 13. Variation of the Fourier parameters during the observed time span, calculated for subsets of 2 d duration each. Left 
panels: at the top, the amplitude Ai is displayed to indicate where Blazhko maxima and minima are located, while lower panels 
show the amplitude ratios -R2I1 ^31) ^41 > and -R51. Right panels: at the top, the phases (pi, (p2, 933 and (pi are plotted, while the 
lower panels show the phase differences ip2i, 1/231, and (/341 in the sine frame. Where error bars arc not visible, they are smaller 
than the symbols. 



Table 3. Overall modulation parameters of V445 Lyr found 
from the fit to the complete data set for the first 6 harmonic 
orders. 



k 


Rk 




Qk 




1 


2.112 


-0.87 


0.357 


0.005779 


2 


3.179 


0.27 


0.521 


0.001900 


3 


3.991 


-0.24 


0.599 


0.000445 


4 


3.806 


-0.08 


0.584 


0.000114 


5 


5.382 


-0.30 


0.687 


0.000037 


6 


4.715 


-0.39 


0.650 


0.000011 



was only 50 d. Some overlap was allowed to guarantee 
frequency resolution of the Blazhko multiplet. 

It was an important result for the CoRoT star 
105288363 that the phasing between the two types of 
modulation was found t o change (see section 5.3 of 
iGuggenberger et al.ll201ll ). Such a phase change is also 
indicated by the variation of AA^ of V445 Lyr (bottom 
panel of Fig. [12} . 



4-2.3 Fourier parameters 

As the so-called Fourier parameters Rki (amplitude ra- 
tio) and ipki (epoch-independent phase difference) are 
considered useful tools to study and compare the prop- 
erties of RR Lyrae stars, we calculated their time vari- 
ability for V445 Lyr. In contrast to the modulation pa- 
rameters discussed in the previous section, they describe 
the pulsation rather than the modulation properties. The 
Fourier parameters are defined as Rki = A^/A-i and 
'Pki ~ '^k — kip\. The data were subdivided into 239 bins 
of 2 d duration, therefore containing about four pulsa- 
tion periods. On such short timescales, the Blazhko effect 
is expected to play only a minor role. The effect of pe- 
riod doubling, however, causes a large error in some bins 
which are more affected by this effect than others. A fit 
including the fundamental pulsation and 10 harmonics 
was calculated for each bin, and the results are displayed 
in Fig. [T2] (amplitude ratios are shown in the left panels, 
phase differences in the right panels). The average values 
of the parameters (derived from the complete data set) 
are as given in Table [D 

Some interesting details are immediately obvious: 
while the amplitude A\ of 0.18 mag is quite small com- 
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Figure 12. Modulation parameters of V445 Lyr versus time. 
Tlic top panel sliows the variation of Ri , which is defined as 
the amplitude ratio between the right and the left modulation 
side peak, and the bottom panel displays the variation of the 
power differenee AA^ . The variation of A points towards a 
shift in the phasing between the two types of modulation. 
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Figure 14. ipi versus Ai diagram of V445 Lyr. In this plot 
it can be seen how the various observed cycles differ from 
each other. Both the phase modulation component and the 
amplitude modulation component change over the observed 
time span. Also, their relative phasing changes from cycle to 
cycle. 



Table 4. Average Fourier parameters of V445 Lyr 



parameter 


value 


parameter 


value 


Ai 


0.184 


VI 


5.84 


R21 


0.268 


^21 


2.34 


R31 


0.096 


V31 


4.95 


Rai 


0.033 


¥'41 


1.08 


R51 


0.013 


¥'51 


3.69 



pared to other RR Lyrae stars, but still in the normal 
range, the amplitude ratios Rki are significantly smaller 
than thos e of non-m odulated stars (for comparison, see 
Fig. 6 in iNemec et a l. 2011). An intriguing feature are 
the sharp upward spikes in the R21 variation which oc- 
cur only during some of the observed Blazhko minima. 
When looking at the phases one notes that, while ipi has 
a smooth periodic variation, the phases Lp2, ipa and Lp4 
show a more or less continuous progression (with excep- 
tions in some cycles), leading to apparent phase jumps in 
the v'fci parameters. 

The stability of the results was tested by using also 
other binnings, and the size of the bins was found to 
play only a minor role. 



4.2.4 Loop diagrams 

A good indicator for the relative contributions of phase 
and amplitude modulation, and for the phasing between 
those two types of modulation, is the Ai versus ipi di- 
agram. The resulting loops for V445 Lyr are plotted in 
Fig. 1141 The direction of motion is indicated with arrows. 
Cycles are defined for this purpose as from one maximum 
of Al to the next, except for the beginning and the end 
of the data set where the additional points are added 
to the adjacent cycles. Cycles 1 to 5 correspond to data 
obtained before the gap in the observations, and cycles 
6 to 8 to data after the gap. Even though the contents 



of Fig. [14] are partly redundant with the upper panels 
of Fig. 1131 this representation allows to better compare 
the observed Blazhko cycles. All observed Blazhko cy- 
cles have a different appearance in this diagram, and the 
contributions of amplitude and phase modulation change 
without notable correlation between each other. 



4.3 O-C diagrams and long-term period change 

In addition to the Fourier analysis which reveals the av- 
erage Blazhko period, and the time-dependent analysis 
which shows the fundamental mode as a function of time, 
we also constructed an O-C diagram, as it can reveal ad- 
ditional details, especially when it comes to long-term 
period changes. The O-C diagram obtained from all max- 
ima in the Kepler V445 Lyr lightcurve is shown in the 
top panel of Fig. 1151 An intriguing feature is the non- 
sinusoidal variation of the O-C values with a few points 
at very low O-C values in some of the cycles. A closer 
inspection of the phase diagrams at the affected times 
reveals that these drops in O-C are happening at the 
epochs with the very unusual distorted light curve shape 
showing double maxima (the double maxima are illus- 
trated in panels b and d of Fig.[T|). The determination of 
the time of the light maximum becomes ambiguous here, 
depending on the maximum that is chosen. The original 
maximum seems to move "to the left" (causing negative 
O-C values) while a bump on the descending branch gets 
stronger and takes the role of the maximum for the next 
Blazhko cycle. This also explains the rather abrupt tran- 
sition from very low to high O-C values. 

We also performed a Fourier analysis on the O-C 
data to check whether the secondary modulation is also 
present in the phase variation. We clearly found the 
Blazhko frequency /s , in this case 0.0184 d"\ as weU as 
the harmonics 2/s and 3/s in the O-C curve, which are 
introduced by the non-sinusoidal variation of the O-C 
curve. The secondary modulation, with a value of fs 
— 0.0064 , was also directly detected in the Fourier 
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Figure 15. O-C diagram of the complete data set of V445 Lyr 
(top panel), and residuals after subtracting a fit with the 
Blazhko frequency, its harmonics and the secondary modula- 
tion frequency (bottom panel) . To guide the eye, a line through 
the means of the residuals in bins of 10 maxima has been plot- 
ted. 

spectrum, as well as the combination peak fs + fs- 
We note that slight differences in the frequency values 
obtained from the O-C diagram compared to those ob- 
tained from the magnitudes do not come unexpectedly, 
as the phase and amplitude modulations are not strictly 
correlated. 

As the quasi-periodic Blazhko modulation is the domi- 
nant signal in the O-C curve, it is necessary to subtract 
a fit including the above mentioned frequencies from 
the O-C curve to be able to identify long-term changes. 
The residuals of this fit are shown in the bottom panel 
of Fig. [15] and reveal quite clearly a long-term period 
change, ft remains unclear, however, whether it is a 
periodic or a continuous linear period change. 



4.4 The analytic modulation approach 

Due to the limitations of the classical Fourier analy- 
sis discussed in Section UJ we applied a new method of 
analysis to the Ke pler data of V445Lyr, that w as re- 
cently described bv lBenko. Szabo fc Papaxdl (1201 ll . here- 
inafter Ell). In this approach, the amplitude and fre- 
quency modulations are treated similarly to the theory 
of electronic signal transmitting, reducing dramatically 
the necessary number of free parameters. This would be 
especially useful for stars like V445 Lyr where a classical 
Fourier analysis yields several hundreds of combination 
peaks. 

As a first step of such an analysis we have to select the fit- 
ting formula, using Table 1 in Bll. The amplitude mod- 
ulation (AM) with the frequency of /b is evident already 
from the shape of the light curve in Fig. [1] Since the 
envelope curve is nearly horizontally symmetric for all 
Blazhko cycles, the AM of /b can be approximated by 
a simple sinusoidal function (see formulae 20 and 21 in 
Bll). 

The high asymmetry of the multiplet peaks' ampli- 
tudes suggests a frequency modulation (FM) as well. Its 



non-sinusoidal nature is clear both from the frequency 
variation function (panel a of Fig. Ilip and the O — C di- 
agram (Fig. I15|) . The combined sinusoidal AM and non- 
sinusoidal FM modulation of /b can be described by the 
formula (49) in Bll, where g' = 1. 

As mentioned above V445 Lyr shows a secondary 
modulation /s as well which is included in the form of an 
AM because of the changing amplitudes of the Blazhko 
cycles. As a first approximation we also assumed this 
modulation as sinusoidal. The linear combination of /b 
and fs shows interaction between the two modulations, 
therefore we have to apply the formula of modulated 
modulation (AM cascade - eq. 27 in Bll). 

The situation of the FM in fs is a bit controversial. 
The existence of an FM seems to be well-established, 
based upon the detection of fs in the Fourier analysis 
of the O-C diagram; however, a combined AM with FM 
in fs does not improve the significance of our numerical 
Levenberg-Marquardt fit. This may be explained on the 
basis of the long cycle length of this modulation and/or 
its weak FM. 

The used best fit model contains a sinusoidal AM and 
a non-sinusoidal FM represented by a three-term Fourier 
sum for the Blazhko modulation and a sinusoidal AM 
for the second modulation. The two modulations are as- 
sumed to be modulating each other. The free parameters 
are the pulsation frequency fo and its harmonics' am- 
plitude and phase up to the 9th order {Ai, A2, ■ ■ ■ , Ag, 
(p-i,(p2, ■ ■ ■ jifig); the modulation frequency fs, the ampli- 
tudes and phases of its AM (obii Vbi) and FM (oai, fm, 
i — 1,2,3), the secondary modulation frequency fs and 
its AM modulation parameters (flsj, ipsi)- They represent 
32 parameters (with the zero point aoo). The model light 
curve shows the global properties of the observed one (see 
Fig, nop , however, the variance of the residual (observed 
minus fitted) is surprisingly high (0.0025 mag). 

The large variance may have different reasons. Some 
of them are method-specific, others are object-specific. 
An important limit of the method is (as mentioned by 
Bll) that it does not describe the migration of the humps 
and bumps caused by the Blazhko effect, a phenomenon 
which is exceptionally strong in V445 Lyr. The situation 
is demonstrated well in the upper panels of Fig. 1161 The 
other problem is that our method assumes regular signals. 
The light curve of V445 Lyr shows, however, irregular be- 
haviour. The loop diagram in Fig.[l4]illustrates the cycle- 
to-cycle variations of the relative strengths of the AM 
and FM components of the modulation. Any static (reg- 
ular) models including the Fourier method face similar 
troubles when using them for such a time-dependent (ir- 
regular) phenomenon. We conclude that therefore, even 
though applying the method lead to a success in obtain- 
ing a reasonable fit with a comparably small number of 
parameters, it is nevertheless not optimal for a compli- 
cated data set like the one on V445 Lyr. 
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Figure 16. The light curve fit of V445 Lyr during Q1-Q4 using the modulated signal approach. The bottom panel shows the 
global 32- parameter fit (red continuous line) with the original data (black dots). The upper panels show a local fit (left) and 
residual light curves (blue crosses) in the interval indicated by the rectangle in the bottom panel. 



5 SPECTROSCOPY, COLOUR 

PHOTOMETRY AND FUNDAMENTAL 
PARAMETERS 

In the framework of ground-based follow up observa- 
tions, three spectra were obtained with the HIRES 
spect rograph at the Keck 10 m telescope in August 
2011 |n emec. Cohen fc Sesaij[201^ ) . Exposure times were 
1200 s, and according to the ephemerides given in sec- 
tion l4. li the spectra were obtained shortly after maximum 
light. Due to its faintness (_Kp=17.4 mag), V445 Lyr is 
not an easy target to observe. From a preliminary anal- 
ysis, heliocentric radial velocities of -392, -390 and -388 
km/s were obtained, as well as a me t allicit y of [Fe/II]=- 
2.0 ± 0.3 dex on the ICarretta et al] l|2009l ) scale, corre- 
sponding to Z=0.0003. 

Based on the metallicity value derived from spec- 
troscopy and with the frequency value of the second over- 
tone mode, we were able to constrain the mass and lumi- 
nosity of V445 Lyr with the help of a theoretical Petersen 
dia gram which is shown in Fig. 1171 Linear pulsation mod- 
els (jSmolec fc Mosk alik 2008) were calculated for a set of 
masses (0.55, 0.65 and 0.75 Mq, which are plotted as 
solid, dashed and dotted lines, respectively) and different 
luminosities (40, 50, 60, 70 Lq, plotted as circles, pluses, 
squares and triangles, respectively). The metallicity val- 



ues necessary to theoretically fit the observed frequencies 
agree very well with the measured metallicity, and the lu- 
minosity and mass of V445 Lyr found from Fig. [17] are 
40-50 Lq and 0.55-0.65 Mq, respectively. 

As the model period ratios depend on the opacity 
tables and the abundance mixture used in the com- 
putations, we tested the stability of ou r results. The 
effect s of different opacit i es an d mixtures (|Asplund et al.l 
I2009I : [g revesse fc Noeld I1993I ). however, were checked 
and found to play only a minor role. In Fig. 1171 the 
models based on O PAL opacities and the mixture of 
lAsplund et all ()2004l ) are shown. V445 Lyr is plotted as 
a black square in the diagram, while the other RRab 
stars for which the presence of the second overtone has 
been reported, are shown as open squares. From the 
wide spread which the stars show in Fig. [17] it is obvious 
that very different parameters are needed to model the 
different stars that show the second overtone. Larger 
masses are neccesary to fit the period ratios of the stars 
with longer periods, while higher metallicities are needed 
to to obtain a model for the shorter periods. Models for 
V350 Lyr and V3 54 Lyr have already be en shown by 
iBenko et all (|2010l ) in their Fig. 6, as well as lNemec et al.l 
( 20111 ) in their Fig. A2. As it was already noted in sec- 
tion [4713] that the stars in which the second overtone is 
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Figure 17. Petersen diagram for V445 Lyr (shown as black 
square) for the secon d overtone, ba sed on models computed 
with the Warsaw code jSmolec fc Mo skalik 2008). Other RRab 
Lyr stars with possible second overtones (as discussed in sec- 
tion |4!T3J are plotted as open squares. 



excited form a very diverse sample, it is not surprising 
that they also cover a wide range in mass and metallicity. 
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Figure 18. Same as FigureHbut for CoRoT 105288363. Note 
that the multiplet structure is more symmetrical than in V445 
Lyr, in the sense that approximately the same number of mul- 
tiplet components appear on both sides of the fundamental 
mode. The peaks belonging to the secondary modulation are 
in similar positions as in V445 Lyr. Again, several additional 
peaks occur which cannot be identified as combination peaks 
of either one of the modulations. This is most likely due to 
the either irregular or yet unresolved changes of the Blazhko 
effect. Unlike in V445 Lyr, no combinations including a 2fg 
term could be found. 



5.1 Ground- based color photometry 

To complement the Kepler data with colour information 
on our target, V445 Lyr was observed from the ground 
in BVRI (of which the RI bands are in Cousins system) 
using telescopes at the Lulin Observatory, including the 
Lulin One-meter Telescope (LOT) and Lulin 0.4m SLT, 
and the Tenagra II Observatory (TNG, with a 0.8m tele- 
scope). The imaging data were reduced with IRAF in a 
standard manner, including bias and dark subtraction, 
and flat-fielding. Phot ometry was obtained fro m the im- 
ages using SExtractor (jBertin fc Arnoutd l 19961) ■ and cal- 
ibrated to the sta ndard m a gnitu des using standard stars 
observation from iLandold (|2009l ). Further details about 
the telescopes, the CCDs and th e reduction of th e imag- 
ing da ta can be found in ISzab6~e t al. (20l3) and iNgeowl 
(|201ll ). who used the same instrumentation for monitor- 
ing the Cep heid V1154 Cyg lo cated within the Kepler 
field-of-view (|Szab6 et al.ll201ll '). 

The observations were performed between 29 March 
and 24 July, 2011 (i.e., during the course of two Blazhko 
cycles) and comprise of about 70 measurements per 
filter with typical uncertainties of about 0.06 mag. As 
the data cover the pulsation cycle well, they allow the 
determination of an average brightness in each colour. 
The following average magnitudes in the standard 
system were obtained for V445 Lyr on the basis of the 
magnitudes of the single measurements: 5=17.80 mag, 
y=17.38 mag, 7?=17.09 mag and 7=16.81 mag. 



6 COMPARISON WITH COROT 105288363 

With the tools developed for the analysis of V445 Lyr, 
we revisited CoRoT 105288363 to apply the same tech- 
niques in a consistent way. Unlike in previous studies, 
the frequencies were kept as free parameters during the 
prewhitening and fitting procedure, in spite of the in- 
crease of computing time. This has the advantage that 
the Blazhko period can be determined not only from 
one measured distance between 2 peaks, but from a 
large number of independently detected peaks. The stan- 
dard deviation of that set of measured values also gives 
a good error estimate. The Blazhko period found by 
this method is 34.6 ± 1.1 d and our solution agrees 
within the error of the prev ious published value of 35.6 d 
(|Guggenberger et a^]|2011^ . 

When the echelle diagram diagnostic was applied 
to CoRoT 105288363, some previously undiscovered fea- 
tures could be unveiled. The first and most important 
finding is a well-resolved secondary modulation which is 
very similar to the one in V445 Lyr, in the sense that 
it has a ratio with the primary modulation of /s//s 
= 2.5 ± 0.27, which is close to the value of /s /fs = 
2.7 ± 0.12 in V445 Lyr. Also, its combination peaks ap- 
pear in similar positions: They are clearly found in the 
higher-frequency side of the harmonics, and also preferen- 
tially on the higher-frequency side of the classical Blazhko 
multiplets. The echelle diagram for CoRoT 105288363 
is shown in Fig. 1181 One has to note that, unlike in 
V445 Lyr, no systematic deviation of the harmonics from 
their expected positions (see section 14.1.2^ is found in 
CoRoT 105288363. Therefore, no tilt of the orders in the 
echelle diagram can be seen. There are, however, some 
irregularities in the pattern with some peaks showing de- 
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Table 5. Comparison of the main characteristics of the two 
stars. Fundamental mode frequencies are very different, as well 
as the Blazhko periods and the periods of the secondary mod- 
ulation. Interstingly, however, the ratio between the primary 
and the secondary modulation are very similar. 
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harmonic order k 





CoRoT 105288363 


V445 Lyr 


/o 


1.76231 


1.94903 


^'o [d] 


0.5674 


0.5131 


Pb [d] 


34.6 


53.1 


/fl[d-l] 


0.0289 


0.0188 


/s[d-i] 


0.0115 


0.0069 


Ps [d] 


86.5 ± 9 


143.3 ± 5.8 


Pb/Ps 


2.5 ± 0.27 
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Pb/Po 


60.9 


103.5 



Figure 19. Same as Figure[5]but for CoRoT 105288363. Here, 
in contrast to V445 Lyr, the components belonging to the 
secondary modulation never reach higher amplitudes than the 
classical multiplet, which might be connected to the fact that 
the changes in the Blazhko effect of CoRoT 10528836 are not 
so dramatic as they are in V445 Lyr. Their amplitude decrease 
is almost as rapid as for the classical components. 
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Figure 20. Fourier spectrum of CoRoT 105288363 after sub- 
traction of 135 frequencies around the fundamental mode and 
its harmonics. A zoom into the region between /g and 2/o is 
shown, clearly revealing the four additional peaks discussed in 
the text. The former positions of the multiplets around /q and 
2/o are marked with shaded boxes. 



leads to the disappearance of all aliases) significantly re- 
duced the noise level of the residuals in the frequency re- 
gi on of interest. The second o vertone that was reported 
bv lGuggenberger" et al.l (|201lh was also found in the new 
analysis (see Fig. 1201 where it is labeled /2). Due to the 
reduced noise, also combinations of /2 with fc/o up to 
fc=5 could be detected, clearly indicating that the fre- 
quency is not due to noise, and not caused by a back- 
ground star. Moreover, three additional frequencies could 
be found (see also Fig. [20}. A peak at 2.3793 d"^ with 
an amplitude of 0.3 mmag is found to be most likely the 
first overtone due to its ratio of /o//i = 0.741 with the 
fundamental mode. Also, a resonance with /o and /2 is 
possible in this case, because 2/i — /2 — /o is only 0.01. An- 
other peak appears at /jv=2.4422 d~^ with an amplitude 
of 0.35 mmag. It cannot be attributed to any overtone. 
This and the fact that it appears between the positions 
of the first and the second overtone makes it quite similar 
to the observations in V445 Lyr (see section |4.1. 5 |l . Also, 
its amplitude is slightly higher than that of the suspected 
first overtone, as is the case for V445 Lyr. The third fre- 
quency found in this range is at /jv, 2=2. 2699 d~^ and has 
an amplitude of 0.3 mmag. It cannot be identified as a 
radial overtone, and it has no counterpart in V445 Lyr. 

Table [5] compares the main characteristics of the two 
stars, including the ratios of the Blazhko modulation, 
while Table [6] compares the additional (overtone) modes. 

In CoRoT 105288363, there is no sign of period 
doubling, and no frequencies could be detected at or near 
the positions of the characteristic half-integer frequencies 
(3/o/2, 5/o/2, etc.). 



viations from the exact position, though not in a system- 
atic manner. The decrease of amplitudes of the secondary 
modulation components with increasing harmonic order 
is shown in Fig. 1191 Those figures should be compared 
to the corresponding ones for V445 Lyr (Figs. [4] and O 
respectively). 

The data of CoRoT 105288363 were then also in- 
spected carefully for signal outside the vicinity of the 
Blazhko multiplets to find evidence for overtones and 
other possible additional modes. The prewhitening of not 
only the classical Blazhko multiplet but also the sec- 
ondary modulation and the additional peaks (which also 



7 SUMMARY AND CONCLUSIONS 

7.1 An unusual star with unusual phenomena 

V445 Lyr is an RRab star with such a strong Blazhko 
modulation that at Blazhko minimum the peak-to-peak 
amplitude decreases down to 0.07 mag compared to ap- 
proximately 1 mag during Blazhko maximum (a differ- 
ence by a factor of 14), leading to the rather low overall 
amplitude Ai of 0.18 mag. 

The light curve around Blazhko minimum shows a strong 
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Table 6. Comparison of the overtone modes and the addi- 
tional frequency of the two stars. 

CoRoT 105288363 V445 Lyr 



/2[d-i] 


2.9856 


3.3307 


A(/2)[mmag] 


0.5 


1.4 


/0//2 


0.590 


0.585 


/i[d-M 


2.3793 


2.6676 


A(/i)[mmag] 


0.3 


2.7 


/0//1 


0.741 


0.731 


/iv[d-i] 


2.442 


2.7719 


^(/iv)[mmag] 


0.4 


3.6 


fo/fN 


0.722 


0.703 



distortion witii a secondary maximum, making it impos- 
sible to even identify V445 Lyr as an RRab pulsator when 
observed only during Blazhko minimum. 
In V445 Lyr the full variety of all the recently discovered 
new features in RR Lyrae stars - period doubling, strong 
irregular changes in the Blazhko effect, a secondary mod- 
ulation, radial overtone pulsation as well as a non-radial 
mode - are combined in one single star. Therefore, it 
serves as an example of how ultra-precise and uninter- 
rupted space photometry can change our view on seem- 
ingly well-known types of stars. The Fourier phases show 
an unusual behaviour which has not been detected be- 
fore in an RR Lyrae star. Also, the distinct spikes in the 
temporal variation of the Fourier amplitude ratios are a 
previously unknown feature which is most likely caused 
by the pronounced double maximum at those phases. 

7.2 Secondary modulation, irregular behaviour, 
and long-term period change 

In V445 Lyr, we find the most extreme variations of the 
Blazhko modulation known so far. This is partly, but not 
fully, explained by the secondary modulation of 143 d 
which we found in the Kepler data. Irregular/chaotic 
changes of the Blazhko modulation and/or even longer 
modulation periods also seem to be present, leading to 
a dense spectrum of peaks around the harmonics of the 
fundamental mode, in which with classical methods up 
to 771 frequencies would be found. 

The amplitudes of the peaks connected to the secondary 
modulation were found to decrease less steeply with har- 
monic order than the components of the classical Blazhko 
multiplets. This interesting feature yet awaits a physical 
explanation. 

A long-term period change is also present, but it could 
not yet unambiguously be determined whether it is a pe- 
riodic or a linear change (or neither). Future Kepler ob- 
servations in the upcoming quarters will certainly reveal 
more about this long-term change. 

7.3 Additional frequencies 

We find at least four additional frequencies not connected 
to the fundamental mode, its harmonics and the Blazhko 
peaks. One of these peaks was interpreted as the sec- 



ond radial overtone, one could possibly be the first radial 
overtone, one was found to be due to period doubling, 
and the fourth was attributed to a non-radial frequency. 
The second overtone is not always present during the ob- 
servations. A strict dependence of the /2 amplitude on 
the Blazhko phase could not be found. Instead, it seems 
to vary rather irregularly. 

Amplitudes and frequencies of all additional peaks 
change notably during the time span of the data. We con- 
sider it possible that fluctuations in the frequency values 
of /i lead to transient resonances which temporarily ex- 
cite the second overtone. 

The additional peaks form numerous combinations with 
fo and the Blazhko multiplets, (including quintuplet 
peaks) and also with the peaks belonging to the sec- 
ondary modulation, indicating that they are all intrin- 
sic to the target. Altogether, 80 peaks were found above 
the significance level at or near the combinations of those 
4 frequencies with the other intrinsic frequencies of the 
target. 

7.4 Spectroscopy, Petersen diagrams, and an 
alternative method of light curve analysis 

Spectroscopy with the Keck telescope revealed a metal- 
licity of [Fe/H]=-2.0 ± 0.3, and Petersen diagrams calcu- 
lated from linear pulsation models point towards a mass 
of 0.55-0.65 A/0 and a luminosity of 40-50 Lq. 
We also applied the new analytic modulation technique 
to the light curve, and found that the best model con- 
tains a sinusoidal AM, a non-sinusoidal FM, as well as 
a sinusoidal AM for the secondary modulation. Due to 
the migration of a strong bump feature and due to the 
irregular/stochastic changes, however, the method faces 
similar troubles as the Fourier analysis. 

7.5 Comparison with another peculiar star 

A revisit of the data on CoRoT 105288363 revealed a 
secondary modulation period with a similar period ratio 
(2.5) with the primary modulation period as V445 Lyr 
(2.7). 

The new analysis of the CoRoT 105288363 data also 
points towards the excitation of more additional modes 
than the previously published second overtone. A non- 
radial mode as well as the first overtone might also be 
excited. 

V445 Lyr also shows a change in the phasing of the two 
types of modulation (amplitude and phase modulation), 
similar to what was observed in CoRoT 105288363. 
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